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Edge magnetoplasmons in a wide armchair graphene ribbon with a weak superlattice 
potential: finite frequency gaps and zero group velocity 
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We show strong effects of a weak and smooth, on the magnetic length, superlattice potential 
upon edge magnetoplasmons (EMPs) at the armchair edge, with a smooth steplike electrostatic 
lateral confining potential, of a wide graphene channel in the v — 2 quantum Hall effect regime. 
The superlattice potential leads to essential enlargement of a number of EMPs, descend from two 
fundamental EMPs in the absence of superlattice. For the wave vector k x within the first Brillouin 
zone, the EMPs show as the regions of acoustical or quasi-acoustical dispersion, with a finite value 
of group velocity, so the regions with frequency gaps, where a group velocity is nullified at some k x . 
We obtain that for k x — > only for two EMPs the frequency tends to zero as for other EMPs it 
obtains finite values. Strong dependence of dispersion relations of the EMPs from the period of the 
superlattice ao and the distance d from a metallic gate is shown; in particular, for typical size of a 
gap, for characteristic value of the frequency and k x at which the group velocity is reduced to zero. 
At the frequency that corresponds to zero group velocity of pertinent fundamental EMP branch the 
response of the system should present a strong resonance. 

PACS numbers: 73.43.Lp,73.22.Pr,68.65.Cd 



I. INTRODUCTION 

Graphene after experimental discovery of its high- 
quality freestanding samples^ has attracted a strong 
attention^ Charge carriers in a single-layer graphene 
possess a gapless, linear spectrum close to the K and 
K' points^— and manifest behavior of chiral massless 
particles with a "light speed" equal to the Fermi ve- 
locity, vp. Graphene shows a lot of unusual effects, 
e.g.: the Klein paradox 2 '-^, i.e., the perfect transmission 
through arbitrarily high and wide barriers upon normal 
incidence (as far as a Dirac-type Hamiltonian is valid), 
a half-integer quantum Hall effect (QHE) 2 -^— , and a 
zitterbewegung ) 2 ' 9 ' 10 i.e., effect induced by a lateral con- 
finement of Dirac fermions. Properties of the latter effect 
are essentially modified by a strong magnetic field. Ex- 
tra Dirac points in the energy spectrum for superlattices 
in graphene have been obtained if the amplitude of peri- 
odic potential is sufficiently large while its period is small 
cnoughJsir— In particular, this leads to new properties 
of the QHE,ii 

Graphene's edges have also been studied 
considerably^^— ' 14 ' 15 in particular, in connection 
with the QHFj 2 ^— ; for some phenomena it matters a 
type, the armchair or zigzag, of edges££~— Edge mag- 
netoplasmons (EMPs) in graphene have been studied 
only recently^ it is shown that in the v = 2 QHE 
regime at the armchair edge, and in the presence of a 
smooth steplike electrostatic lateral confining potential, 
the chirality, spectrum, spatial structure, and number of 
the fundamental EMPs depend strongly on the position 
of the Fermi level Ep ■ 

In the case (i) of Ref. [l6J, when Ep intesects (see 
Fig. 1, cf. with Fig. 1(a) of Ref. 16) four degener- 
ate states of the zero LL at one location and two de- 
generate states of this LL at a different location, two 



fundamental EMPs are present: counterpropagating and 
with essential spatial overlap. This is in contrast with 
EMPs in conventional two-dimensional electron systems 
(2DES) which give only one fundamental EMP at the 
v = 2 QHE regime, with negligible spin-splitting; for 
conventional 2DES different types of EMPs have been 
studied theoretically^— and experimentally 24 ~ 2 6 Above 
two counterpropagating EMPs can be on resonance if a 
strong coupling of the EMPs holds at the ends of the 
segment L™ % < L x , where L x is the length of graphene 
channel^ 

In present study for the case qualitatively outlined in 
Fig. 1 (i.e., it is the case (i) of Ref. [la) we explore 
theoretically effect of a weak and smooth superlattice 
potential V s (x) — V s cos(Gx) with G = 2n/a , upon 
EMPs. Here, in agreement with speculations o£i& that 
a strong Bragg coupling is possible due to a weak super- 
lattice along the edge (with period L e x m , if L x /L™ 1 > 1), 
we show that V s {x) can have a strong effect on two fun- 
damental EMPs leading to manifestation of resonance 
effects; in particular, referred to in the abstract. Present 
EMPs in graphene with the superlattice are very differ- 
ent from the EMPs treated previously for conventional 
2DES with a superlattice^ 

In Sec. II A we obtain the wave functions and the 
spectra of LLs in an infinitely large graphene flake in 
the presence of a perpendicular magnetic field and of a 
smooth electrostatic confining potential, along the y di- 
rection, as without V s {x) so in its presence. In Sec. II B 
we study the combined effect of a smooth, step-like elec- 
trostatic confining potential and of armchair graphene 
edges, at y = ±L y /2, and of the superlattice potential 
V s (x) on the local Hall conductivity in the v = 2 QHE 
regime. In Sec. Ill we obtain strong renormalization of 
the EMPs in graphene by a weak superlattice potential. 
We make concluding remarks in Sec. IV. 
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Figure 1: (Color online) Energy spectrum of the n — 0, ±1 
LLs as a function of the quantum number j/o, at the right half 
of symmetric graphene channel with armchair edges and the 
smooth electrostatic potential, Eq. ((9J, for the Fermi level 
Ef — Vb/2. Spatially separated edge states are created at y" 
and at yf. (marked by an upward arrow) as the branches of 
the n = LL cross Ef- The v — 2 QHE is manifested in dc 
magnetotransport. 



II. GRAPHENE CHANNEL AND LOCAL HALL 
CONDUCTIVITY 



A. Effect of a smooth potential and of a weak 
periodic potential on the LLs 

We consider a long and a wide flat graphene flake of 
length L x and width L yi with armchair edges at y = 
±L y /2 in the presence of a perpendicular magnetic field 
B = Bz, of a smooth confining potential V y = V(y) 
along the y direction of electrostatic origin, and of a one- 
dimensional (ID) periodic potential V s (x) along the x 
direction. We assume that V s {x) = V s cos(Gx) is a weak 
ID modulation potential of period aa. For definiteness 
we assume that the potential V y is symmetric. 

If it is not otherwise stated, we consider solutions with 
energy and wave vector close to the K point; we present 
pertinent results for the energies and wave vectors close 
to the K' point (valley) as well. In the nearest-neighbor, 
tight-binding model the one-electron Dirac Hamiltonian, 
for massless electrons, is H = Ho + tV s (x) where Ho = 
Vpo ■ p + tV y , with 1 the 2x2 unit matrix. Explicitly 
Ho is given by (e > 0) 



H = v F 



P. i 



V v /vf Px 

- ip y - eBy 



- ip y - eBy 

Vy/VF 



(1) 



where p x and p y are components of the momentum oper- 
ator p and vp ~ 10 6 m/s the Fermi velocity. The vector 
potential is taken in the Landau gauge, A = (—By, 0, 0). 



1. Landau levels for a smooth potential V(y) 

First, we present properties of the LLs in the absence 
of periodic potential, when only a smooth potential and 
armchair termination are assumed^. The equation (Ho — 
E)ip = admits solutions of the form 



V> (0) (r) 



c Hv)/yJL x 



•<»' " (£$) • < 2 » 



where the components &A(y) and &s(y) correspond to 
the two sublattices and the coefficients A and B satisfy 
the relation \A\ 2 + \B\ 2 = 1; r = {x,y}. Introducing 
the magnetic length £o = (H/eB) 



1/2 



</<> 



— ^o^xai the 



variable £ = (y—yo)/^Oi an< i assuming that V y is a smooth 
function of y, with a characteristic scale Ay 3> (-Oi it 



follows 16 at r 
energy E"' 



ioa/hvp ^C 1 for the n 
V(yo) and 



LL that the 



Al 






(i) = V« 
(l/r)[l-«(l 



2M/4 e -?72 



) 1/2 ]£°, 



(3) 



where a — dV/dt; |^ = o, and K = +(— ) corresponds to the 
K (K') valley, here B% = A°_ ta 1 and A° + = B°_ ^ 
r/2<l. 

Finally, without periodic potential for any n — 
0,±1,±2, ... LL and yo not too close to the graphene 
lattice termination at y = ±L y /2 (see Fig. 1), the eigen- 



values E, 



(0),K 



= E, 



(0), 



can be written as 



EiT = S9 n(n)^ 
to 






V(y ), 



0,±1,... (4) 



where the sign function sgn(n) = 1 and —1 for n > and 
n < 0, respectively. Notice that each n^0 LL is twice 
degenerate with respect to the valley quantum number 
K. However, the latter index is kept in the left hand 
side of Eq. (UJ) as for yo(k xa ) close to the armchair edge 

these eigenvalues, and especially strongly E l' K , become 
dependent on k; notice, for these conditions k can not be 
related only to one valley 6-8 . Accordingly, for n ^ LL 
and yo n °t too close to the graphene lattice termination, 
see Fig. 1, the eigenvalues ([4J are four times degenerate. 
Wave functions pertinent to the eigenvalues Eq. ((4]) we 
will denote also as ?/r ^ (r). 



2. Effect of potential V(y) and of a weak smooth periodic 
potential V s (x) on the LLs 

Assuming hv F /£o > hv g (k xa )G > V s /2 and G£ < 1 
where the group velocity v g (k xa ) = hr^-dE^l /dk xa , 
now we will study how a weak periodic potential V s (x) 
modifies energies and wave functions of the LLs for a 
smooth potential V(y); in particular, for conditions per- 
tinent to Fig. 1. We calculate the eigenvalues and 
eigenfunctions corresponding to the Hamiltonian H = 



Wo + tV s (x) using the perturbation theory 28 . Because 
further there will be important only the eigenstates lo- 
calized, along y, near a right edge of the channel, we will 
write formulas assuming, in particular, such eigenstates. 
Similar with Refill, we can neglect by a small "nonres- 
onance" contributions, np ^ n a . Then keeping only the 
"resonance" contributions, np = n a , and the terms of 
the first order over V s , e.g., for the eigenfunctions of the 
n = LL we obtain 



r 0<kxa (r)=4± K a 



V„ 



2hv g (k xa ) G 



^(0),« 



.^(0),« 



Further, the eigenvalues are well approximated by the 



(0), 



zero-order terms, i.e., E* k sa E y n £ , where, for j/o not 

too close to the graphene lattice termination, E t K are 
given by Eq. (@}. Indeed, the first order corrections are 
exactly nullified and the second order ones are very small, 
e.g.: ~ (V S /2V ) 2 < 1 for the model potential Eq. ©. 



can be well approximated as the valley index. A strong 
splitting between the electron, n — +, and the hole, k = 
— , branches of the n = LL 2 ^^, due to hybridization of 
the valley states take place nearby the armchair edge, at 
\L y /2 — yo\ < la- The eigenvalues of the n = LL for k — 
+(—) increase (decrease) with increasing yo. However, 
for the n > 1 LLs the k = ± branches at the armchair 
edge have a small splitting, due to hybridization of the 
valley states, as their eigenvalues increase with increasing 
2/o ; these branches are attributed to the electron band. 
Notice, the electron, (n = 0, k = +), LL stems from the 
conduction band and the hole, (n = 0, k = — ), LL arises 
from the valence band. 

We now consider the situations depicted in Fig. 1 for 
a wide symmetric armchair graphene ribbon L y > 2y r 3> 
Ay ^> £ . For definiteness the smooth lateral potential 
is assumed as follows 

V(y) = (V /2)\2+my-y r )/Ay) + $((y+y r )/Ayj\, (9) 



B. Local Hall conductivity in the v = 2 QHE regime 

Extending magnetotransport formulas for the local 
Hall conductivity a yx (r) of a standard 2DES in the chan- 



29-31 



nel, in the presence of a smooth lateral potentia 

obtain, for linear responses and in strong magnetic fields, 

a yx (r) in the form 3 - 2 



o-yx(r) = n(r)e/B, 



(0) 



where the local electron density n(r) is smooth, on the 
characteristic scale £o, as along y, mainly monotonic, so 
along x, with a weak periodic modulation. It is given by 



n(r) 



an 



f aK (a>K\16(r - r)\aK), 



(7) 



with a = {n,k xa }; a yx {y) = -a xy (y). So far in Eqs. 
©-([7]) only the electrons from the conduction band LLs 
are assumed. However, if the valence band LLs can es- 
sentially contribute to a yx (r), as for Fig. 1, the local hole 
density p(r) will contribute to the right hand part of Eq. 
© by changing n(r) on [n(r) — p(r)}. Then, for condi- 
tions relevant to Fig. 1, when only (n = 0, k = ±) LLs 
can essentially contribute to a yx (r) or a diagonal com- 
ponent of the local conductivity tensor, equation (JB]) is 
rewritten as 



&yx(x,y) = 



2e 2 L x 



/OO 
dyo[fo, Vo , K - S K ,-} 
-oo 



X(l>5,k.j*)\\1>0,k. a (*)), 



(8) 



where f„ 



is the Fermi function, the two-component 



^o,fc XQ ( r )) is given by Eq. 



column spinor wave function 

([5|); the factor 2 accounts for spin degeneracy. In Eqs. 
Q),© k — ± is understood as the pseudospin quantum 
number—; e.g., for y > only at (L y /2 — yo)/^o ^> 1 it 



where $(x) is the probability integral. In Fig. 1 we have 
L y = 18Ay, V = hv F /V2£ , y r = 5Ay, and Ay = 10l . 
When the Fermi level Ep is between the bottoms of the 
n = and n = 1 LLs, at yo — 0, and the condition Vq 3> 
2UbT holds, the occupation of the n > 1 LLs is negligible; 
the same holds for the n = LL in the regions of yo 
that are well above Ep, see Fig. 1. In addition to the 
smoothness of the potential Eq.®, we assume armchair 
edges of the graphene sheet at y = ±L y /2, which cause 
the bending of the LLsr^^ and L y /2 — y r > Ay. For 
conditions of Fig. 1 and qualitatively similar, the dc 
magnetotransport measurements will manifest the v = 2 
QHE. 

Further, for conditions qualitatively similar with those 
of Fig. 1, in agreement with speculations of RefJ£ we 
will show that a weak periodic potential V s (x) can have 
a strong effect on two fundamental EMPs leading to man- 
ifestation of the resonance effects. First, for convenience 
of a reader, we will present expressions for the local Hall 
conductivity, obtained in Ref. 1 ^, that are pertinent to 
V.(x)->0. 



1. Effect of a smooth potential, an armchair edge and a 

weak periodic potential on local Hall conductivity in the 

v — 2 QHE regime 

Now, for different regions of the graphene channel we 
will present expressions as for, obtained in ReLi&, the un- 
perturbed local Hall conductivity, a yx (y) , so for the main 
contribution induced by a finite V s (x), <r yx (x 7 y) ex V s . 

The superscript in a yx (x, y) indicates that it is of the 
first order over V s . Correspondingly, this contribution in 
Eq. ((5J) stems from the first order contributions to the 
wave function Eq. ([5]); to calculate a yx (x, y) we will need 
also assume that Vs/ksT <§C 1. Point out that a yx {y) 
is given by the right hand side of Eq. © if to substi- 



tute tpQ k (r) \ , calculated by taking into account the 
first order corrections, by the zero order wave function 

In case (i), for y > and {yf — yo)/£o 3> 1, from Eqs. 
©-©, ©-© h follows^ 



4°J(y) 



2e 2 



tanh 



YM1 - YM 

2k B T 



(10) 



where V(y) is so smooth on the scale of £q that 
£ dV(y^)/dy <C k B T; the factor 4 accounts for spin 
and pseudospin degeneracy. Introducing the charac- 
teristic length £t = £o(k B T£o/hv^), this condition of 
smoothness can be rewritten as £q <C £t, where u" = 
£\hr x dV {%£) J dy is the group velocity at the edge y". 
Notice, for conditions of Fig. 1 it follows that v^/vp — 

(£ /V2irAy) <g; 1, due to £ /Ay < 1. Point out, Eq. 
(fTU)) shows that at y — y" the Hall conductivity changes 
its sign from the electron type of charge carriers to the 
hole one. 

For (£ T /Ay) 2 < 1, Eq. {jTUJ) can be rewritten asM 



4°J(v) 



h 



tanh 



Vr -y 

2£ T 



From Eq. $TB) it follows!^ 

dv { yx (y) 



dy 



4e^ 
ft 



Mi 



cosh 



2£ T 



(11) 



(12) 



Further, for L y /2 > y > L y /2 — 5£q pertinent numerical 
results ^ 14 ' 15 for v = 2 we approximate by the same 
analitical expression for [n(y) — p{y)] as in Ref. Qi|, That 
gives 



A°)(oA = 



2e 2 



< y J(y) = — 



dyn 



-(y-yoV/t 2 o 



[fo 



.yo . 



-1], (13) 



where it is used that Eq^ Vo> - is a sharply decreasing func- 
tion at yo ~ vf such that the Fermi function in Eq. (p~3|) 
is very fastly growing at yo « yf on a scale £d ^ £o- Here 
appears a new characteristic scale £d = /csT^g//i|tig| as 
for a change of yo on £j, at yo ~ yf , the value of E . 



will change on fc^T; point out that 
in Fig. 1. 

From Eq. ([13]) it follows 1 ^ 



yo, 



<C vf is implicit 



da y °J(y)/dy=(2e 2 /hV^£ )e 



-(y~y d r ) 2 /e 2 o 



(14) 



by changing the derivatives over y to those over j/o and 
integrating by parts. 

In a similar manner, for case (i) and y > 0, we obtain 

from Eqs. ©-©, © that <r$ (r) = cr$ (y) cos (Gz) , 
where 



*&?(!/) 



v; 



ft, fcgT 



cosh 



2k B Tto ( v _ 



^h\v d q \ 



{y-y d r) 



V-Vr 

2£ T 



(15) 



here it is taken into account that tr < 0. Due to 
fesTfo/ftl^gl <C 1, the amplitude of the second term in 



the square brackets of Eq. ([15)) is less than of the first 
one. In addition, Eq. ([T5|) shows that a small parame- 
ter V s /k B T <C 1 warrants that |<7yx (#, y) <C \(T VX i(y)\', in 
particular, at y ~ y". 



III. STRONG RENORMALIZATION OF THE 

EMPS IN GRAPHENE BY A WEAK 

SUPERLATTICE POTENTIAL 

Now we will study effect of the superlattice potential, 
V s (x), on the fundamental EMPs for case (i), i.e., for 
conditions like in Fig. 1; dissipation is neglected. For 
V s (x) = 0, our fundamental EMPs will coincide with ones 
obtained in Ref. [l6J. Similar with Ref. [lg, we expect that 
the charge excitation due to EMPs at the right part of 
channel will be strongly localized at y" ( p ru (t, r)) and y d 
( p rd (t, r)). Then the components of the current density 
j(w,r) in the low-frequency limit uj <C Vf/£q are given, 
cf. with^ as 

j x ( W ,r) = -[*$>&) +*$(T)]E y ( U ,T) 

+ v u g p™{u,v)+v d g p rd {u,r), (16) 

j y (w,r) = [a y °J(y)+a y 1 ](r)}E x (u J ,r), (17) 

where we have suppressed the factor exp(— icut) common 
to all terms in Eqs. ([j"6[) and (ITTl) . From Eqs. ([T6)) - 
(|17[) . Poisson's equation, the linearized continuity equa- 
tion and using that for EMPs E(w,r) = — V<^(w,r)), we 
obtain 

-iuOffa r) + p rd (u, r)) + [v«d x p™( U) r) 

+^9 x/ / rf (u;,r)] + [9 lff W(r)]S v »,( W ,r) 

-ldy<? y Q J (y) + dytfj (r)]dM", «0 = 0, (18) 

where c^ = d/dx. Point out, in Eq. (ITS)) the coefficients 
are invariant for translations along x on any distance 
integral of the period clq. Then in Eq. (|18p we assume 
that 



p ru {oj,r)= J2pt U (u>k x ,y)e ik y*, 

i 
p rd (u;,r) = J2 p r A",k*,y)e ik * )x , 



e=-i 
l 



?( w , r ) = ^2 <Pi( UJ ,k x ,y)e 



ik^x 



(19) 



where £ = — 1, 0, 1, A4 = &x + 27r£/ao, and fc^ = k x ■ 

In Eq. (TT9")) . for a metallic gate at a distance (i from the 
2DES (e.g., it can be a heavily doped Si separated from 



the graphene sheet by a SiC>2 layer of thickness d = 300 
nm), is given by 



<P£{u,k x ,y) = - / dy'R g (\y-y'\,k x e) ;d) 

e J —OO 

x[p r t u {w,k x ,y')+pl d {u,k x ,y')], (20) 

where R g (---) is given by 
R g (\y-y'\,k^;d) = go(l*ff ) lll/ -y'l) 

- ^o(l4"V(y-y') 2 + 4d 2 ) ; (21) 

where i^o(2 ; ) is the modified Bessel function. Without 
of a metallic gate, d — ¥ oo, the dielectric constant e is 
spatially homogeneous if not stated otherwise. 

Multiplying Eq. (1181) by i exp(—ik x x) and then inte- 
grating over x, J Q x dx, we obtain 

(u - k x v^)p r u {w, k x , y) + (u- k x v d )p r d (w, k x ,y) 
—cri l J(y)[9y<P-i(u,k x ,y)-d y tpi{u,k x ,y)] 



six unknown functions: p r e u {oj,k x ) and p^ d (u>,k x ), where 
t = -1, 0, 1. They read, with y^ u = yf. - y™ > 0, 

[w - w^oCftw; d)K"(w, fc x ) - 2c h k x R g (y du , k x ;d) 



xpg (w,fc x ) + c h k x 



k B T 



Ri(y? u ,k x -V;d) 



x p'J^u, fe x ) + <(yf , fcW; d)^V, fc x )l = 0, (26) 






[w - w^ (fc x ; d)}pl d (uj, k x ) + c h k x R g (y du , k x ; d) 
/ V v u \ r 

xp?i(w, fcx) + R' g {yt u , k x ^;d)p\ u {uj, k x )] = 0, (27) 
xR'Jy? u ,k x ;d)p r d (u>,k x )=0, (28) 



[^»4x(y)]^^o(w, fc x , y) + -[d y <T${y)\ 
k x -\^, k x ,y) + *£V(«, *., »)1 - 0. (22) '" : ^ (4- X) ; <0]p3(«> M + onk^ 



Further, multiplying Eq. (fT8| by iexp(— ik x x) and 
then integrating over x, L x dx, we obtain 

(w - fci" 1 )^)^^, fc xjtf ) + (« - fc^v^plKw, fc„ y) 
— ^ 1 J(y)d y <po(u,k x ,y) + [d v a^J(y)]ki'^tp-i(u;,k x ,y) 



xR g (y du , k x -V;d) P r _ u 1 (u J , k x ) - c h k x ~V 



Vs 

knT 



X£t 2\^\ I R 9 {y r> kx; ^^^ kx) = °' (29) 



oo 



w-«« (fc«;d)]pr(«,fcx)-2c fc *W 



+ 2^[9 !/ crW(j/)](po(w,fc :E , 2 /) = 0. 



(23) xi^f^VKV^+c^ (^«r) 



In addition, multiplying Eq. (J18I) by iexp(— i/c x x) and 
then integrating over x, J„ x dx, we obtain 



X Jl^.fe; (*)/#*(«, fc B )=0, 



W -^ (feW;d)]^ d (u;,fe x ) + cafci 1) 



(30) 



no 



+ -k x [d y (T$(y)}(po(uJ,k X7 y) = 0. 



(24) 



xR g (y d \k x V;d)pr(LU,k x )-c h k x V 



i)( V * 



knT 



For yf - ifi > ^ T , from Eqs. ([2"2 > (P ]l it follows that 
p r e u (to , k Xj jj) and p rd (y> 1 k x ,y) can be well approximated 
(cf. Ref. M) as 

^(w.fe,!/) = U^cosh^^-^)] -1 ^™^,^), 
^(w.fe.i/) = (l/V^4)e- fe -^ )2/£ V^(c,fc ;c ).(25) 

In addition, we can neglect by overlap between 
p ru {uj,k x ,y) and p rd (uj,k x ,y) in Eqs. p2]>-(|23]>. Then, 
by integration of Eqs. (I2"2"j) - (|2"41 over y within separate 
regions around y" and y d , we obtain by straightforward 
calculations six coupled linear homogeneous equations for 



eT 2\^\ ) W"' kx] d)/9 ™ (w ' kx) = °' 



(31) 



where, by using the same notations as in Ref. Il6l . we 
have Ch = 4e 2 //ie, 

u>+ (k x ,d) = k x Vg +2chk x a p (k x ;d), (32) 



uj_ (k x ,d) = -kx\Vg\ - c h k x a m (k x ;d), (33) 
with the matrix elements 



a p (k x ;d) = — 



OO /*OC 



t&rcft i? 5 (£T|^ — f|, &#; d) 



2/j./o\ ' 



oo •/ — oo 



cosli (x/2)cosli (£/2) 



(34) 
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: 2 ~~ 2 ~ " 


-1 

-2 
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Figure 2: (Color online) The dispersion relations u>(k x ,d — ► 
oo), without the gate, of six EMPs calculated from Eqs. (|26|l - 
500cm" 1 , Vs/kaT = 0.3, y, d " = 104, Ay = 



?/to 



2, »: 



4 x 10 6 cm/s, 



r,; - —3 x 10' cm/s, e = 2. Panel (a) presents the dispersion 



(131)1 for 27r/a 

104, B = 9T, T = 77K, 

d , „ ln 7 

U 9 - 

relations within the first Brillouin zone, n/ao > k x > — 7r/oo, 
by the curves 1 (solid), 2 (dashed), 3 and 4 (dotted), 5 and 6 
(dash-dotted). Panel (b) presents a zoom of the anticrossing 
for the branches 2 and 4, at k x w —165 cm -1 and w ~ — 1.36x 
10 11 s~\ with the gap » 0.98 x 10 9 s _1 ; here the EMPs 2 and 
4 have a zero value of group velocity for pertinent k x ■ Panel 
(c) presents a zoom of the branches 1 and 2 at k x « and 
a; ~ 0; here a finite gap is absent. Panel (d) presents a zoom 
of the anticrossing for the branches 3 and 5. 



Q'm.y.fcx] d) 



— oo J — oo 



dxdt 

P x 2 +t 2 



R g (£ \x-t\,k x ;d). (35) 



Point out, if the graphene is located between two dielec- 
tric media with dielectric constants: ei for the halfspace 
below of the graphene, and £2 for the dielectric layer be- 
tween the graphene and a metallic gate, - then in all 
above expressions e = (ei + £2)/2. 
In addition, in Eqs. pH) - (j3Tj) 



R' g (y*-y?,k x ;d) = \k x \[K 1 (\k x \y d r u ) 

.da 



!), 



V%?") 2 + 4d 2 



K^k^ {y^Y + Ad?) 



(36) 



Notice that in the long- wavelength limit, k x £r <C 1, and 
for large d, such that the effect of gate, oc exp(— 2|fc x |d) -C 
1, can be neglected, from Eqs. (|2"Tj) . ([32l) -(|36 ]) we obtain 
that: R g (\y-y'\,k x ;d) w ln(2/\k x {y- y')\) -7, where 7 
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Figure 3: (Color online) The dispersion relations ui(k x ,d = 
3000nm) of six EMPs calculated from Eqs. pB > l[3"T ]l for 
2ir/ao — 500 cm -1 and other parameters of Fig. 2 except 
d = 3000nm. Panel (a) presents the dispersion relations 
within the first Brillouin zone by the curves 1 (solid), 2 
(dashed), 3 and 4 (dotted), 5 and 6 (dash-dotted). Panel 
(b) presents a zoom of the anticrossing for the branches 2 and 
4, at k x « —189 cm -1 and u> ~ — 1.17 x 10 11 sT 1 , with the gap 
~ 0.90 x 10 9 s^ 1 . Panel (c) presents a zoom of the branches 
1 and 2 at k x « and u) ~ 0. Panel (d) presents a zoom of 
the anticrossing for the branches 3 and 5. 



is the Euler constant, a p (k x ;d) « [ln(l/|fca;|#r) — 0.145], 
a m (k x ;d) w [ln(l/\k x \e ) + 3/4], and R' g (y? u , k x] d) w 

l/y*». 

For ^ = 2 and case (i) , in Fig. 2 we plot the dispersion 
relations u>(k x , d —> 00) of the EMPs calculated from Eqs. 



Ay 



U for 2Tr/a = 500 cm" 1 , V s /k B T = 0.3, 
104, B = 9T, T = 77K, ^ T /£ = 2, ^ 



</ /;. 



= 4x 
10 B cm/s, »J = -3 x 10 7 cm/s, e = 2, and £ w 8.5 
nm. Fig. 2(a) presents the dispersion relations of these 
EMPs within the first Brillouin zone by the curves 1-6. 
Notice, Fig. 2(a) shows that the branches 3 and 4 have 
to w ±2.0 x 10 n s _1 at k x w 0. Fig. 2(b) presents a 
zoom of the anticrossing for the branches 2 and 4, at 
fcj, rj —165 cm" 1 and ui ~ — 1.36 x 10 11 s~ 1 , with the gap 
sw 0.98 x 10 9 s" 1 : here the EMPs 2 and 4 have a zero 
value of group velocity for pertinent k x . A panel for the 
anticrossing of the branches 1 and 3 it follows from the 
Fig. 2(b) by changing k x on — k x and co on — w. Fig. 
2(c) presents a zoom of the branches 1 and 2 at k. x w 
and w « 0; here a finite gap is absent. Panel (d) presents 
a zoom of the anticrossing for the branches 3 and 5 at 
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Figure 4: (Color online) The dispersion relations ui{k x ,d = 
300nm) of six EMPs calculated from Eqs. p6 )l -(f3l ]l for 
2iv/ao = 500 cm -1 and other parameters of Fig. 2 except 
d = 300nm. Panel (a) presents the dispersion relations within 
the first Brillouin zone by the curves 1 (solid), 2 (dashed), 3 
and 4 (dotted), 5 and 6 (dash-dotted). Panel (b) presents 
a zoom of the anticrossing for the branches 2 and 4, at 
k x « -222cm" 1 and w « -0.90 x 10 u s -1 , with the gap 
~ 0.73 x 10 9 s" 1 . Panel (c) presents a zoom of the branches 
1 and 2 at k x « and u) « 0. 



k x « -159.4cm" 1 and w « 2.611 x 10 n s -1 , with the 
gap « 1.37 x 10 6 s" 1 . 

In Fig. 2 and below, it is used that the exact dispersion 
relation ui{k x ,d) of any EMP mode can be presented in 
the form periodic in the reciprocal space, i.e, ui{k Xl d) = 
Lo(k x ± 2ir/ao,d), and continuous across the borders of 
the Brillouin zone, w(7r/ao — 0, d) = w(7r/ao + 0, d). The 
latter, in particular, does not allow an infinite group ve- 
locity for the EMP. Point out that the dispersion curves 
1 (solid), 2 (dashed), 3 and 4 (dotted) have correct peri- 
odic and continuous form in the reciprocal space, k Xl and 
both qualitatively and quantitavely well describe disper- 
sion of pertinent EMP modes in graphene with the su- 
perlattice: 1 and 2 are the main fundamental EMPs, 3 
and 4 are the first excited fundamental EMPs. However, 
an approximate dispersion curves 5 and 6 (dash-dotted) 
qualitatively correctly represent pertinent exact depen- 
dencies only nearby the anticrossings of 5 with 3 and of 
6 with 4. So the curves 5 and 6 are shown only within 
a small part of the first Brillouin zone in Fig. 2(a). In 
addition, as the second order contributions over the pe- 
riodic potential are neglected (as well as an additional 
contributions in Eqs. (19) with the £ = ±2) the Fig. 
2(d) gives only rough approximation for this anticrossing 
and, in particular, for its gap. 

In Fig. 3 we plot the dispersion relations u>(k x ,d = 
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Figure 5: (Color online) The dispersion relations ui(k x ,d — ► 
oo), without the gate, of six EMPs calculated from Eqs. (|26|> - 
(|31[) for 2-k/cio = 5000 cm" 1 , the rest of parameters are the 
same as in Fig. 2. Panel (a) presents the dispersion relations 
within the first Brillouin zone, ir/ao > k x > —ir/ao, by the 
curves 1 (solid), 2 (dashed), 3 and 4 (dotted), 5 and 6 (dash- 
dotted). Panel (b) presents a zoom of the anticrossing for 



the branches 2 and 4, at k x 



-1860 cm and lj ; 



-1.15 x 



10 s , with the gap » 9.0 x 10 s. Panel (c) presents a 
zoom of the branches 1 and 2 at k x fa and uj ~ 0; here a finite 
gap is absent. Panel (d) presents a zoom of the anticrossing 
for the branches 3 and 5. 



3000nm) of the EMPs calculated from Eqs. ([2"6 > (|3l ]l for 
d = 3000nm, the rest of parameters are the same as in 
Fig. 2. Fig. 3(a) shows that the branches 3 and 4 have 
lj « ±1.9 x lO 1 ^" 1 at k x sw 0. Fig. 3(b) presents a 
zoom of the anticrossing for the branches 2 and 4, at 
k x fa —189cm" 1 and u ~ —1.17 x 10 11 s" 1 , with the gap 



0.90 x 10 a 



here the EMPs 2 and 4 have a zero 



value of group velocity for pertinent k x . A panel for the 
anticrossing of the branches 1 and 3 it follows from the 
Fig. 3(b) by changing k x on — k x and lo on — u. Fig. 3(c) 
presents a zoom of the branches 1 and 2 at k x sa and 
uj sa 0; here a finite gap is absent. Fig. 3(d) presents 
a zoom of the anticrossing for the branches 3 and 5 at 
k x fa -121cm" 1 and u ~ 2.34 x 10 11 s _1 , with the gap 
~5xl0 5 s _1 . 

In Fig. 4 we plot the dispersion relations ui(k x ,d = 
300nm) of the EMPs calculated from Eqs. ([26 ]) - ([3l] ) for 
d = 300nm, the rest of parameters are the same as in 
Fig. 2. Fig. 4(a) shows that the branches 3 and 4 have 
u ~ ±1.6 x 10 11 s _1 at k x sa 0. Fig. 4(b) presents a 
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Figure 6: (Color online) The dispersion relations uj(k x ,d = 
3000nm) of six EMPs calculated from Eqs. p5 | -l[gT )l for the 
gate at d = 3000nm; as 2ir/ao — 5000 cm" 1 so the other 
parameters are the same as in Fig. 5. Panel (a) presents the 
dispersion relations within the first Brillouin zone. Panel (b) 
presents a zoom of the anticrossing for the branches 2 and 4, 
at k x ps —1920 cm" 1 and w ~ —1.13 x 10 12 s _1 , with the gap 
with the gap ps 8.8 x 10 9 s _1 . Panel (c) presents a zoom of 
the branches 1 and 2 at k x w and u) w 0; here a finite gap is 
absent. Panel (d) presents a zoom of the anticrossing for the 
branches 3 and 5. 



zoom of the anticrossing for the branches 2 and 4, at 



and u> 
-l. 



-0.90 x lO^s" 



with the 



k x ps -222 cm- 
gap ps 7.3 x 10 8 s _1 : here the EMPs 2 and 4 have a zero 
value of group velocity for pertinent k x . A panel for the 
anticrossing of the branches 1 and 3 it follows from the 
Fig. 4(b) by changing k x on —k x and uj on —uj. Fig. 
4(c) presents a zoom of the branches 1 and 2 at k x ps 
and uj ps 0; here a finite gap is absent. The anticrossing 
for the branches 3 and 5 holds at k x ps —56.7 cm" 1 and 
wps 1.79 X lO 1 ^" 1 . 

In Fig. 5 we plot the dispersion relations uj(k x , d — > oo) 
of the EMPs calculated from Eqs. flU-flS]} for 27r/a = 
5000cm" 1 by the curves 1-6. Fig. 5(a) presents the dis- 
persion relations of these EMPs within the first Brillouin 
zone. Notice, Fig. 5(a) shows that the branches 3 and 4 
have uj ps ±1.787 x 10 12 s _1 at k x ps 0. Fig. 5(b) presents 
a zoom of the anticrossing for the branches 2 and 4, at 



1860 cm- 



and 



1.15 x 10 12 s"\ with the 



k ES 

h, x IV 

gap ps 9.0 x 10 9 s _1 : here the EMPs 2 and 4 have a zero 
value of group velocity for pertinent k x . A panel for the 
anticrossing of the branches 1 and 3 it follows from the 
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Figure 7: (Color online) The dispersion relations u)(k x ,d = 
300nm) of six EMPs calculated from Eqs. ([26])- ([31]) for 
2ir/ao = 5000 cm" 1 and other parameters of Fig. 5 except 
d = 300nm. Panel (a) presents the dispersion relations within 
the first Brillouin zone. Panel (b) presents a zoom of the an- 
ticrossing for the branches 2 and 4, at k x ~ —2220 cm" 1 and 
lj ps -8.95 x 10 11 s" 1 , with the gap ps 7.3 x 10 9 s" 1 . Panel (c) 
presents a zoom of the branches 1 and 2 at k x ~ and w ~ 0; 
here a finite gap is absent. 



Fig. 5(b) by changing k x on — k x and uj on — uj. Fig. 
5(c) presents a zoom of the branches 1 and 2 at k. x ps 
and ui (=a 0; here a finite gap is absent. Panel (d) presents 
a zoom of the anticrossing for the branches 3 and 5 at 
k x ps —1112cm" 1 and u ps 2.16 x 10 12 s~ 1 , with the gap 
~5xl0 6 s -1 . 

In Fig. 6 we plot the dispersion relations u>(k x ,d = 
3000nm) of the EMPs calculated from Eqs. pg| ) -(f5T ]) for 
2-7r/a = 5000 cm" 1 . Fig. 6(a) presents the dispersion 
relations of these EMPs within the first Brillouin zone by 
the curves 1-6. Notice, Fig. 6(a) shows that the branches 
3 and 4 have uj ps ±1.783 x 10 12 s" 1 at k x ps 0. Fig. 6(b) 
presents a zoom of the anticrossing for the branches 2 and 
4, at k x ps -1920cm" 1 and w ps -1.13 x 10 12 s _1 , with 
the gap ps 8.8 x 10 9 s" 1 : here the EMPs 2 and 4 have a 
zero value of group velocity for pertinent k x . A panel for 
the anticrossing of the branches 1 and 3 it follows from 
the Fig. 6(b) by changing k x on —k x and ui on — u>. Fig. 
6(c) presents a zoom of the branches 1 and 2 at k x ps 
and u; ps 0; here a finite gap is absent. Fig. 6(d) presents 
a zoom of the anticrossing for the branches 3 and 5 at 
k x ps -1081cm" 1 and wps 2.143 x 10 12 s"\ with the gap 
~3.7x lO^" 1 . 

In Fig. 7 we plot the dispersion relations to(k x ,d = 
300nm) of the EMPs calculated from Eqs. ([26 ]) - ([3l] ) for 
d = 300?im, the rest of parameters are the same as in 
Fig. 5. Fig. 7(a) shows that the branches 3 and 4 have 
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Figure 8: (Color online) The dispersion relations ui{k x ,d — ► 
oo), without the gate, of six EMPs calculated from Eqs. (|26|l - 
(J3TJ for 2-R/ao — 50000 cm -1 , the rest of parameters are the 
same as in Fig. 2. Panel (a) presents the dispersion relations 
within the first Brillouin zone, ir/ao > k x > — ir/ao, by the 
curves 1 (solid), 2 (dashed), 3 and 4 (dotted), 5 and 6 (dash- 
dotted). Panel (b) presents a zoom of the anticrossing for 
the branches 2 and 4, at k x Rs — 22150 cm -1 and to ~ —8.6 x 
10 12 s _1 , with the gap « 6.9 x 10 10 s _1 . Panel (c) presents a 
zoom of the branches 1 and 2 at k x Rs and u Rs 0; a finite 
gap is absent. 



u ~ dzi.ou x iir~s - at K x « 0. Fig. 7(b) presents 
a zoom of the anticrossing for the branches 2 and 4, at 
k x « -2220cm -1 and to « -8.95 x 10 n s -1 , with the 
gap « 7.3 x 10 9 s _1 : here the EMPs 2 and 4 have a zero 
value of group velocity for pertinent k x . A panel for the 
anticrossing of the branches 1 and 3 it follows from the 
Fig. 7(b) by changing k x on —k x and to on — lo. Fig. 7(c) 
presents a zoom of the branches 1 and 2 at k x « and 
w sw 0; here a finite gap is absent. 

In Fig. 8 we plot the dispersion relations Lo{k x , d —} oo) 
of the EMPs calculated from Eqs. P5 |) -([3l | for 27r/a = 
50000 cm -1 . Fig. 8(a) presents the dispersion relations of 
these EMPs within the first Brillouin zone by the curves 
1-6. Notice, Fig. 8(a) shows that the branches 3 and 4 
have to ~ ±1.42 x 10 13 s -1 at k x w 0. Fig. 8(b) presents 
a zoom of the anticrossing for the branches 2 and 4, at 
k x m -22150 cm -1 and to ~ -8.6 x 10 12 s"\ with the 
gap « 6.9 x 10 10 s -1 : here the EMPs 2 and 4 have a zero 
value of group velocity for pertinent k x . A panel for the 
anticrossing of the branches 1 and 3 it follows from the 
Fig. 5(b) by changing k x on — k x and lo on — to. Fig. 8(c) 
presents a zoom of the branches 1 and 2 at k x « and 
lo k, 0; here a finite gap is absent. 

In Fig. 9 we plot the dispersion relations Lo(k x ,d = 
3000nm) of the EMPs calculated from Eqs. pJt - fBT]) 
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Figure 9: (Color online) The dispersion relations u>(kx,d = 
3000nm) of the EMPs calculated from Eqs. p6 > l(3T ll for pa- 
rameters of Fig. 8, except d = 3000nm. Panel (a) presents 
the dispersion relations within the first Brillouin zone. Panel 
(b) presents a zoom of the anticrossing for the branches 2 and 
4, with the gap ~ 6.9 x 10 10 s -1 . Panel (c) presents a zoom 
of the branches 1 and 2 at k x (a and u) f» 0; a finite gap is 
absent. 



for d = 3000nm; in addition, 2ir/a,Q — 50000 cm -1 and 
other parameters coincide with those of Fig. 8. Fig. 
9(a) presents the dispersion relations of these EMPs 
within the first Brillouin zone by the curves 1-6. No- 
tice, Fig. 9(a) shows that the branches 3 and 4 have 
lo w ±1.42 x 10 13 s -1 at k x rj 0. Fig. 9(b) presents 
a zoom of the anticrossing for the branches 2 and 4, at 
k x fa -22150 cm -1 and to sa -8.6 x 10 12 s -1 , with the 
gap « 6.9 x 10 10 s -1 : notice, the parameters of this anti- 
crossing are very close to the ones of similar anticrossing 
in Fig. 8(b). A panel for the anticrossing of the branches 
1 and 3 it follows from the Fig. 9(b) by changing k x on 
— k x and lo on —to. Fig. 9(c) presents a zoom of the 
branches 1 and 2 at k x f=a and lo « 0; here a finite gap 
is absent. 

In Fig. 10 we plot the dispersion relations Lo(k Xl d = 
300nm) of the EMPs calculated from Eqs. p6 ]) -([3l ]l for 
d = 300nm; in addition, 2ir/ao = 50000 cm -1 and other 
parameters coincide with those of Fig. 8. Fig. 10(a) 
presents the dispersion relations of these EMPs within 
the first Brillouin zone by pertinent curves, 1-6. No- 
tice, Fig. 10(a) shows that the branches 3 and 4 have 
lo rj ±1.41 x 10 13 s -1 at k x rj 0. Fig. 10(b) presents 
a zoom of the anticrossing for the branches 2 and 4, at 
k x rj -22900 cm" 1 and lo rj -8.2 x 10 12 s -1 , with the gap 



6.7x10 



10, 



Notice, the parameters of this anticross- 



ing are close to the ones in Fig. 8(b), in particular, in 
Fig. 10(b) the gap is only 3% smaller than in Fig. 8(b). 
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Figure 10: (Color online) The dispersion relations ui{k x ,d — 
300nm) of the EMPs calculated from Eqs. J26)-([3TJ for pa- 
rameters of Fig. 8, except d — 300nm. Panel (a) presents the 
dispersion relations within the first Brillouin zone. Panel (b) 
presents a zoom of the anticrossing for the branches 2 and 4, 
at k x w -22900 cm -1 and to « -8.2 x 10 12 s _1 , with the gap 
~ 6.7 x 10 °s _1 . Panel (c) presents a zoom of the branches 
1 and 2 at k x « and 10 ~ 0; a finite gap is absent. 



A panel for the anticrossing of the branches 1 and 3 it 
follows from the Fig. 10(b) by changing k x on —k x and 
u on — u). Fig. 10(c) presents a zoom of the branches 
1 and 2 at k x ~ and uj rj 0; here a finite gap is ab- 
sent. Notice, Fig. 10(a) shows that the anticrossing of 
the curves 3 and 5 takes place at k x ~ 1730 cm" 1 and 
u « 1.45 x 10 13 s -1 . 

It is natural to call in Figs. 2-10: the fundamental 
EMPs 1, 2 as the main fundamental EMPs, and the 
fundamental EMPs 3, 4 as the first excited fundamen- 
tal EMPs. Figs. 2-10 show that for case (i), outlined 
in Fig.l, a strong Bragg coupling is possible due to a 
weak superlattice along the edge, with the period ao, if 
L x /ao ^$> 1. In particular, for frequencies in the THz 
range: cf. Figs. 8-10. We expect that for the frequency 
that corresponds to zero group velocity of pertinent main 
fundamental EMP branch, or pertinent first excited fun- 
damental EMP branch, and its vicinity the response of 
the system will have a strong resonance. 



IV. 



CONCLUDING REMARKS 



At the edge of a wide armchair graphene ribbon in 
the v = 2 QHE regime and with a smooth monotonic 
electrostatic potential, we investigated the appearance of 
EMPs that show zero group velocity, for characteristic 
frequencies, and finite frequency gaps due to effect of 



a weak superlattice potential. The superlattice poten- 
tial from two original fundamental EMPs (present here 
without superlattice^), due to a strong Bragg alike cou- 
pling of them, gives correctly within as a whole first Bril- 
louin zone two main fundamental EMPs (branches 1 and 
2 on Figs. 2-10) and two first excited fundamental EMPs 
(branches 3 and 4 on Figs. 2-10). In addition, for the 
wave vector within the center of the first Brillouin zone, 
i.e., k x — > 0, only the frequencies of the main fundamen- 
tal EMPs, 1 and 2, tend to zero as the frequencies, e.g., 
of the first excited fundamental EMPs, 3 and 4, tend to 
finite values. As at the frequency that corresponds to 
zero group velocity of the main fundamental EMPs so at 
one for the first excited fundamental EMPs the response 
of the system should have a strong resonance, e.g., in the 
THz range; see Figs. 8-10. 

Next we list and discuss the approximations used. 
Point out that in Fig. 1 (as well as in Fig. 1 o£i&) it is im- 
plicit that (v g (k xa )/vF) 2 <C 1 for any shown yo — tnk xa . 
Then extra Dirac points in the energy spectrumiir— due 
to present smooth, weak superlattices will not appear as 
here it follows that TwfG ^> V s /2, which is opposite to 
the key condition for extra Dirac pointsi 11 ' 12 Here for 
the EMPs in the v = 2 QHE regime dissipation is ne- 
glected, which is well justified as here the EMP damping 
can be related only with inelastic scattering processes 
within narrow temperature belts, of width fc^T, of each 
edge state that are much weaker than scattering processes 
due to a static disorder^ The latter makes a dominant 
contribution to the transport scattering time in a 2DES 
of grapheneii£i22 for B = 0. Needless to say that for a 
more accurate account of the EMPs studied here, dissi- 
pation must be included in the treatment. We have ne- 
glected by nonlocal effects that usually have minor effect 
on fundamental EMPs 23 . We emphasize that our study 
of the fundamental EMPs for the armchair termination 
of a graphene channel cannot be directly extended to 
zigzag termination as some important properties of the 
wave functions and the energy levels are different than 
those of the armchair termination, cfi 2 -^— ^. We rele- 
gate the study of EMPs along zigzag edges to a future 
work. 

It is used a simple analytical model of a smooth, lateral 
confining potential Eq. ([§]), however, our main results 
are robust to modifications of its form and parameters 
if the qualitative conditions of Fig. 1 are realized in a 
graphene channel in the v = 2 QHE regime. In Figs. 
2-10 it is used that the exact dispersion relation u(k x , d) 
of any EMP mode can be presented in the form periodic 
in the reciprocal space, i.e, u>(k x ,d) = u>{k x ± 27r/ao,d), 
and continuous across the borders of the Brillouin zone, 
ui(ir/ao — 0, d) = ujijr/ao + 0,d). The latter, in partic- 
ular, docs not allow an infinite group velocity for the 
EMP. Point out that the dispersion curves of first four 
EMP modes, 1 - 4 , have correct periodic and continuous 
form in the reciprocal space, k x , and both qualitatively 
and quantitavely well describe dispersion of these EMP 
modes in graphene with the superlattice. However, an 



11 



approximate dispersion curves 5 and 6 (dash-dotted) in 
Figs. 2-10 qualitatively correctly represent pertinent ex- 
act dependencies only nearby the anticrossings of 5 with 
3 and of 6 with 4. So the curves 5 and 6 are shown only 
within a small part of the first Brillouin zone in Figs. 
2(a)-10(a). In addition, as the second order contribu- 
tions over the periodic potential are neglected (as well as 
an additional contributions in Eqs. (TK)|) with the t = ±2) 
the Figs. 2(d), 3(d), 5(d), 6(d) give only rough approx- 
imation for this anticrossing and, in particular, for its 
gap. 
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